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Abstract 

Motivated by recent data on B ^ ppK decay, we study various charmless three-body baryonic B 
decay modes, including Apir, S'^pvr, ppvr, ppK^ , in a factorization approach. These modes have rates 
of order 10^^. There are two mechanisms for the production of baryon pairs: current-produced 
and transition. The behavior of decay spectra from these baryon production mechanisms can be 
understood by using QCD counting rules. Predictions on rates and decay spectra can be checked 
in the near future. 

PACS numbers: 13.25.Hw, 14.40.Nd 
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I. INTRODUCTION 



The Belle collaboration recently reported the observation of B ppK decay, the first 
ever charmless baryonic B decay mode, giving B = {A.Stol ± 0-5) X 10-6 [ig. Three-body 
baryonic decay in 6 ^ c transitions have been observed previously, following a suggestion 
by Dunietz It is interesting to compare the charmless case to the charmful one and also 
to charmless two-body modes such as B^ pp, which has B < 1.2 x lO^^ |Q. 

It has been pointed out that reduced energy release (e.g. by a fast recoiling meson) would 
favor the generation of baryon pair and thus three-body baryonic modes could be enhanced 
over two-body rates One of the signatures would be (baryon pair) threshold enhance- 
ment in the three-body baryonic modes. In our previous study of B^ D*^pn 0, we 
assumed factorization and obtained up to 60% of experimental rate from the vector current 
contribution. The decay spectrum exhibits threshold enhancement. The same threshold 
enhancement effect was predicted for the charmless ppn mode, giving B ~ 10"^ 0. It is 
interesting that the newly observed ppK mode shows such a threshold enhancement [|l|. 
With this encouragement we extend our study to charmless modes such as Apn, Spyr, ppir' , 
ppK~ and ppK^ . These modes are interesting not just by their (possibly) large rates, but 
also by their accessibility. Some of these modes are studied in a recent work that utilizes 
a factorization and pole model approach. 

In Sec. II, we extend the factorization approach to the charmless case, where one now has 
two mechanisms for baryon pair production. In Sec. Ill, we discuss baryonic form factors 
and their associated Quantum Chromo-Dynamics (QCD) counting rules In Sec. IV, the 
formulation is applied to the above mentioned charmless modes. The threshold enhancement 
phenomenon is found to be closely related to the QCD counting rules. Discussion and 
conclusion are given in Sec. V, while some useful formulas are collected in an appendix. 



B 




FIG. 1: (a) The current-produced [J) and (b) transition (T) diagrams for B ppK decay. 
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II. FACTORIZATION 



In this section we extend the factorization approach used in Refs. p, |^ to charmless 
decay modes. Under the factorization assumption, the three-body baryonic B decay matrix 
element is separated into either a current-produced baryon pair {J') part together with a B 
to recoil meson transition part, or a i? to baryon pair transition (T) part together with a 
current-produced recoil meson part. As an example, the current-produced and transition 
diagrams for B~ ppK~ decay are depicted in Fig. 

In charmless decay modes, we need to use the effective Hamiltonian consisting of operators 
and Wilson coefficients, which is standard and can be found, for example, in Refs. ||T^, IT 



In this work, we concentrate on the dominant terms. The factorization formula for the decay 
process B ^ XY, where (X, Y) = {h, BB ) with h being a light meson and BB a baryon 
pair or vice versa is given by 

G ( 

M{B->XY) = --^iKfeKJctil^w)^-^ ® (Mfe)y-A + a2(MM)y-A ® 



-VtbVtl [ag ^{q'q')v~A ® {qb)v-A + ^4 ^{qq')v-A ® {q'b)v-A 

q' q' 

+«5 ^iq'q')v+A ® iqb)v~A - 2a6 ^iqq')s+p ® iq'b)s~p 

q' q' 

+ ^ag ^ Cg'iq q')v~A ® {qb)v-A + . . . ] | , 



(1) 

q' ^ 

where Oi ® O2 stands for {X\Oi\0) {Y\02\B) . The coefficients a^ are defined in terms of 
the effective Wilson coefficients as aj=odd = c"^^ + c^^i/Nc and aj=cven = c^^ + c^^i/^c- 



TABLE I: The coefficients Oj for & — > s [6 — > d] from [10|. Values for og — og are in units of 10 . 





Nc = 2 




iVc = 3 


Nc = 00 


ai 


0.99 [0.99] 




1.05 [1.05] 


1.17 [1.17] 


02 


0.22 [0.22] 




0.02 [0.02] 


-0.37 [-0.37] 


as 


-4.5 - 23 i [-2 - 20 i] 


72.7 


-0.3i [73 + 0.3 


227 + 45 i [223 + 41 i] 


04 


-349.5 - 113.5 i [-338.5 - 101.5 i] 


-387.3 - 


12H [-375.7- 108.3 i] 


-463- 136 i [-450 - 122 i] 


05 


-166 - 23 i [-164- 20 z] 


-66- 


-0.3i [-66 + 0.3 i] 


134 + 45 i [130 + 4H] 


ae 


-533 - 113.5 i [-523 - 101.5 i] 


-555.3 - 


12H [-544.7- 108.3 i] 


-600- 136 i [-588 - 122 i] 


ag 


-86.8 -2.7i [-86.6 -2.5i] 


-92.6 - 


-2.7i [-92.4- 2. 5«] 


-104.3 - 2.7i [-104.1 - 2.5i] 
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We stress that cf^ are renormalization scale and scheme independent, as vertex and penguin 
corrections are included [|T^. Their values are given in Table |. In this work we use the 
iVc = 3 case, while the iY^ = 2, cx) cases are shown to indicate non-factorizable effects. 

III. FORM FACTORS AND QCD COUNTING RULES 

In this section, we first discuss the meson form factors used in this work. We then turn 
to discuss baryonic form factors, especially the implication of QCD counting rules. 

A. Meson Form Factors 



The decay constant fh of the pseudoscalar meson h is defined as 

(%,)|g7^(l-75)g'|0) = ^^p^ 



(2) 



These parameters and quark masses are taken from Ref. ||TT . 
We also need 0^ 0^ form factors defined as follows: 



(hip,) \qY{l-j,)b\B{pB)) 



{PB + PhT - ^^^^-^ iPB - PhT 



{Pb - Phf 
ml -ml 



(Pb - PhY 

We use the so-called MS form factors, which take the following form [1^ : 



{pB - PhY Fi-\t) . (3) 



:pB^h 




;i - t/M^) [1 - t/M^ + 0-2 tyM^] 





l-ait/M^ + a2tyM^ ' 
where My = 5.42 (5.32) GeV for h = K{7c). Other parameters are given in Table |I[ 



(4) 
(5) 



TABLE IL Relevant parameters for the B ^ K,7t transition form factors of Eqs. (H) and dq). 





pB^K 




pB-.. 


pB^TT 

^0 


F,y>^{o) 


0.36 


0.36 


0.29 


0.29 




0.43 


0.70 


0.48 


0.76 


0-2 




0.27 




0.28 
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B. Baryon Form Factors 



Factorization introduces two types of matrix elements containing the baryon pair: 
(BB'|l^(y4)|0) involving vector (V) or axial vector (A) current-produced baryon pair, and 
(B'B'IV {A)\B) involving the B — > BB' transition. 

For the current-produced matrix elements, we have 

me + J 

= ^(pb) + ^2)7, + — (p^, -pb)"| ^^(PbO , (6) 
I m-B + ' ^ ' ) 

(BBV.IO) = ^(pb) (^A(t)7. + ^^^,^2 (pb+Pb'),|75^(Pb'), (7) 
where Fi 2 are the induced vector form factors, qa the axial form factor, and the induced 
pseudoscalar form factor. We have used Gordon decomposition to obtain the second line of 
Eq. (^. Note that t = {pb +Pb')^ ^ "^bb' nothing but the BB pair mass. 

According to QCD counting rules [^, both the vector form factor Fi and the axial form 
factor qa, supplemented with the leading logs, behave as in the t ^ 00 limit, since 
we need two hard gluons to distribute large momentum transfer. F2 and Ha behave as 
acquiring an extra 1/t due to helicity flip. In the electromagnetic current case, the 
asymptotic form has been confirmed by many experimental measurements of the nucleon 
magnetic (Sachs) form factor Gm = Fi + F2, over a wide range of momentum transfers in 
the space-like region. The asymptotic behavior for also seems to hold in the time-like 
region, as reported by the Fermilab E760 experiment |]13[ for 8.9 GeV^ < t < 13 GeV^. 
Another Fermilab experiment, E835, has recently reported for momentum transfers 

up to ~ 14.4 GeV^. An empirical fit of \G^m\ = Ct^'^[ln{t/Ql)]'^, is in agreement with the 
QCD counting rule prediction. 

The current induced form factors Fi + F2 and qa can be related by means of the SU(3) 

decomposition form factors Fv,a, Dv,a and Sv,a (with Sv,a appearing only in the non- 

traceless current shown in Table |IT|. It is well known that Fy and Dy can be 

expressed by the nucleon magnetic form factors G^", 

1 3 
Fv = + ^GIj , Dy = --GIj. (8) 

As the first term Fi + F2 in Eq. (|^) can be related to nucleon magnetic Sachs form factor 
Gm, similarly the second term F2 can be related to {Ge — Gm) /[t/ ("^b + ^b')^ ~ 1]) where 
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TABLE III: Relations of baryon form factors Fi + F2 and qa with the nucleon magnetic form 
factors GK'I^. 



BB 


V, A 


SU(3) 


Fi + F2 


00) 


pn 


{ud)v,A 


(F + D)v,A 




'-J f 'f 1 /TIL 

3 '-^M ^J'M 


Ap 


( — \ 








SOp 


(sn)y,A 




7I i^M + 2 Gil) 




pp 




{F + D + S)v,A 




3 "-^Af '-'Af 


pp 


idd)v,A 


Sv,A 


Sv 


1 /-~fP O/^'T' 

3 "-^A// ^'-"A^ 


pp 


{ss)v,A 


{D-F + S)v,A 







pp 


{uu + dd + ss)v,A 


{2D + 3S)v,A 


3 {-G\j + 'S'y) 




pp 


(cu uu + eddd + es ss)v,a 


{F + D/3)v,A 




^A/ 



G^; is the nucleon electric Sachs form factor. Since we do not have enough data on time-like 
nucleon Ge, we concentrate on the Fi + F2 term as we did in Ref. p. We may in fact gain 
information on G^; by reversing our present analysis on these three-body baryonic B decays 
in the future when more data become available. 

The nucleon magnetic form factors are fitted to available data in Ref. by 



i=l 



V — 



i=l 



0/ J 



(9) 



where 7 = 2.148, xi = 420.96 GeV^ X2 



-10485.50 GeV^ 



X3 



106390.97 GeV^ X4 



-433916.61 GeV^o, = 613780.15 GeV^^^ yi = 236.69 GeV^ y2 = -579.51 GeV^ and 
Aq = 0.3 GeV. They satisfy QCD counting rules and describe time-like electromagnetic 
data such as e+e" A^A^ suitably well. We have real and positive (negative) time-like G^/"'' 



|15| , |T6|. It is interesting to note the alternating signs of the Xi and yi parameters, and that 
only two terms are needed to describe the neutron magnetic form factor [0. 

The time-like form factors related to Sy, Fa, Da, Sa are not yet measured. It is noted 
in Ref. ^ that the asymptotic behavior of baryon form factors studied in the 80s may be 
useful. Their asymptotic behavior ast ^ 00 can be described by two form factors depending 
on the reacting quark having parallel or anti-parallel spin with respect to baryon spin [jl^ . 
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By expressing these two form factors in terms of 



as t 



oo, one has 



Da 



Fa 



Sv 



Sa 




(10) 



Since these relations only hold for large t, it implies relations on the leading terms of these 
form factors. In general more terms may be needed. In analogy to the neutron magnetic 
form case, we express these form factors up to the second term 



The asymptotic relations of Eq. ([T0| ) imply si = xi — 2yi, fi = 2xi/3+?/i/2, di = xi — 3?/i/2 
and §1 = —Xi/3 + 2 yi. 

The coefficients of the second terms are undetermined due to the lack of data. However, 
we can use the axial vector [g^ = Fa + Da) contribution in —>■ D*~pn decay to constrain 
/2 and d2. The part of the branching fraction B{B^ D*~pn) = (14.5l|^ ± 2.7) x 10"^ § 
arising from the vector current has been calculated to give By ~ 7 x 10^'^ [§. We find for 
/2 + c?2 = —2110 GeV^, the branching fraction coming from the axial current Ba{B^ — *• 
D*~pn) ~ 12.7 X 10~^, and the sum By + Ba is within the measurement range. Had we 
used the asymptotic form of Eq. ( |IOD for in the whole time-like region, we would obtain 
Ba ~ 1-0 X 10^'^, which is too small. 

In this work we take 



and S2 = S2 = for simplicity, since there is no data to constrain these yet. It is interesting to 
note that the asymptotic relations give vanishing results for {pP\{ss)v,a\0) (~ Dy^A — Fy^A + 
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(11) 



f^ + d2 = z = -2110 GeV^ 



(12) 



Sv,a) , as one can see from Table |T| and Eq. ([T0|) . This can be understood as OZI suppression. 
We still have a vanishing {pp\{ss)a\0) for smaller t, if we take d2 = f2 = z/2. Since there is 
no point to advocate a large ss form factor in this work and this choice is preferred by the 
OZI rule, we therefore use d2 = f2 = z/2 throughout. For the vector case, we have vanishing 
{pp\{ss)v\0) if we use the asymptotic relation in the whole time-like region. On the other 
hand, if we use Eq. (^) for Sy, we may have a small but non- vanishing (ss)y form factor for 
small t. This may be related to the pole effect in the VMD view point [Q. Furthermore, 
we find that other OZI suppressed current-produced matrix elements, such as {nn\{ss)v,A\0), 
{^+^-\{dd)v,A\0), (S-S+|(MM)y,A|0), {E-E-^\{uu)v,a\0) and (S°H°|((irf)v,A|0), have the same 
SU(3) decomposition as the {pp\{ss)v,a\0) one. They therefore do not provide any further 
constraint. 

We will also encounter (BB \qq'\0), which can be related to the vector matrix element by 



the equation of motion. 



(BB'igViO) = i^^±^(BB'|r,|0) 



niq — rriqi 

rriB — Tn^' ,,,,,, , . 

= — t w t^ PB' • 13 

mq — rriqi 

This gives safe chiral limit in the B 7^ B' case. For example, in (Ap|Ms|0), we have 
(mA — mp)/{ms — my) ~ 1. If B = B', we encounter (me — m^)/{mq — rriq). As hinted 
from the Ap case, our ansatz is to take this factor as the number Uq of the corresponding 
constituent quark in B. For example, we take {pp\dd\0) ~ Fiuv, while {pp\ss\0) ~ as 
suggested by the OZI rule. 

For Ha, we follow Ref. p| to control the behavior of pseudoscalar form factors in the 
chiral limit by using 

\2 



hA{t) = --^ 2^ 9A{t) , 14 

t — mi 



''GB 

where stands for the corresponding Goldstone boson mass. Thus, in the chiral limit. 



(BB'|g75g'|0) = (^^+Pb')^ (BB'|A,|0) 

mq + mq' 



9A{t)u{pBh5v{p^') (15) 



m,q + m„/ t — 



stays finite, otherwise, we will be facing a large enhancement factor (rriB + m^/)/ {mq + m^/) 
in the above equation as we turn off hA{t). 
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We now turn to the transition form factors. In general, the matrix element of 5 — BB 
transition can be defined as 

(BB'lg (75) h\B) =iu{pB) [^v{5) h + ^a{5) ^75 + ^p(5) 75 + ^5(5)] viPB')^ (16) 

where Ph = Pb ~ Pb — Pb', and the form factors J-'v = Ta-o = = -^P5 = by parity 
invariance. According to QCD counting rules, for t = m?—, 00, we need three hard 

BB 

gluons to distribute the large momentum transfer released from the b q transition. An 
additional gluon kicks the spectator quark in the B meson such that it becomes energetic 
in the final baryon pair. Thus, as t — ^ 00, we have: 

1 1 

That J?-p55 have one more power of 1/t than J^a,V5 is due to helicity flip. This can be easily 
seen by taking \Bg) ~ ^75^710). Without any chirality flip due to quark masses, we only have 
J-'a,v5 and the above counting rule holds, while with additional chirality flip, more effectively 
from the b quark mass, we can also have J-'p.s5 but with additional power of 1/t. 

In this work we will need the transition matrix elements {pp\{ub)s,p\B^) and 
{pp\{db)s,p\B^), which consist of eight form factors in total. It is useful to restrict these 



even if only by some asymptotic relations. By following a similar path to Ref. [17], the 
chiral conserving parts J-'a,v5 can be expressed by two form factors depending on the inter- 
acting quark having parallel or anti-parallel spin with the proton spin. The chiral flipping 
parts J^p^s5 can be expressed by one form factor with the spin of interacting quark parallel to 
the proton's. The spin anti-parallel part is absent since it corresponds to an octet-decuplet 
instead of an octet-octet baryon pair final state. The asymptotic forms (as m^, t — » 00) are 

{pp\{ub)s\B~) = iu{pp) [Fa ^hlb + Fp-/^] v{pp), 
{pp\{ub)p\B') = iu{pp) [Fv5 i)h + Fp] v{pp), 
{pp\{db)s\B^) = iu{pp 



^(11Fa + 9Fv5) i>hl^-\Fpi^ 



{pp\{db)p\B^) = iu{p 



V) 



^(QF^ + llFys) h-\Fp 



v{Pp), 

v{Pp). (18) 



We need only three form factors. For simplicity, we use 

Ca.V5 t7 



i^AV5 = ^, FP = ^, (19) 



and Eq. (p!8D in the whole time-like region. 



IV. CHARMLESS BARYONIC B DECAYS 



We now apply the results of the previous sections to charmless B — » ApvT"*", S^pTr"*", ppK , 
and —>■ ppTT~ , ppK~ decays. These modes are of interest not just because of possibly 
large rates, but also by accessibility in detection. 

Let T denote the part of the decay amplitude M. that involves the B — > BB transition 
matrix element (BB |\/(y4)|i?), and JT" denote the part that involves current-produced BB' 
matrix element (BB |V(/1)|0). We have, 

M (A(S0)p7r+) = J (A(S0)j)7r+) , 

Ai {pph) = J (pph) + T {pph) , 

J{ppW) = J {ppK-) , (20) 

where h = K~'^ ^ vr"''^. We note that A(E°)p7r"'" modes only have current-produced contri- 
butions J . On the other hand, the ppTT~ mode is dominated by T{ppTT^) contributions, as 
we will see later. It can be used to extract baryonic transition form factors which can be 
applied to ppK, ppK^ modes via Eq. (pISl). Furthermore, the ppK^ and the ppK~ modes 
have identical current-produced matrix elements, as one can easily show by replacing the 
spectator quark in the B^ K~ transition. 



A. W 



Ap7T+, 



By using Eq. (|I]) and equations of motion we have. 



M(A(E°)p7r 
Gp 



= J{A{J:')pn+) 

7i+\uY{l - 75)615°) ((KfeK>i - Vt,VZa,){BsP\sj^{l - 75)«|0) 



+ 



nis -rriu + m„ 



(21) 



rUb - ruu 

As stated before, only the current-produced part (JT") contributes, hence it is similar to the 
B^ K^*^^-K^ mode, where one only has transition while, analogous to Ap and 

the K^*^~ is produced by the current |TT . 



The chiral limit of the vector term is protected by baryon mass differences, while that of 
the axial term is protected by Eq. (|T3|). Since the contribution from the vector current {V) 
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TABLE IV: B{B^ Ap7r+[S'^j)7r+]) (in units of 10 ^) as decomposed into the vector {By) and 
axial vector {Ba) current contributions. The form factor inputs are explained in the text. 





iVc = 3 


Nc 


= CO 


Bv Ba 


Bv Ba 


Bv 


Ba 


C^j;^ 0.13 [0.51] 0.07 [0.16] 


0.14 [0.56] 0.08 [0.17] 


0.18 [0.66] 


0.10 [0.21] 


Fa, Da — 0.31 [0.27] 


— 0.35 [0.30] 




0.42 [0.36] 




3.5 4 
m( A, z») p ( GeV) 



FIG. 2: dB/dnij^op (upper three lines) and dB/dm\p (lower three lines) plots for 72 = ^2 = z/2. 
Dashed, solid and dot-dashed lines are for Nc = 2, 3, oo, respectively. 

does not interfere with that from the axial current {A), the branching fraction is a simple 
sum of the two, i.e. B = By + Ba- 

Taking 03 (or 7) = 54.8° from a recent analysis [l^, we give in Table |I^ the branching 
fractions for A(S°)p7r+. The first row is obtained by extending the asymptotic 

relations to the whole time-like region and information from the nucleon magnetic form 
factors G^j^P for qa, while we apply f2 = d2 = z/2 in the second row. Note that we do 
not need Sa,v- Our result for the Apn^ mode in the first line is consistent with that of 
Ref. [p. We concentrate on the Nc = 3 case, while Nc = 2, 00 cases are given to indicate 
possible non-factorizable effects. Since B A(S0)p7r+ are penguin-dominated processes, 
the branching fractions are dominated by the ag term of Eq. (pT]). One can verify this by 
comparing different Nc cases in Tables and |. Since Nc dependence is weak in this term, 
we do not expect large non-factorizable contributions. 

The axial contribution to A(S°)p7r"^ mode in the second row is about four (two) times 
larger than that in the first row. We find i3(A(S°)p7r+) ~ [0.5 (0.9)] x 10^^, giving a larger 
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. 4 



45^ 



90° 
<t)3 



135° 180' 



FIG. 3: 03-dependence of the branching fractions of i? Apn^ (lower three hnes) and 
Ti^pTT^ (upper three hnes) for 72 = ^2 = z/2. Notation the same as Fig. ^. 

rate for the S°p7r+ mode. These rates are similar to those obtained in [^f. 

We show in Fig. ^jthe dB/dmp^p and the dB/dmj^op for the f2 = d2 = z/2 case, which gives 
larger rates. One clearly sees threshold enhancement, which can be seen as a consequence 
of the need for large-t suppression of the baryon form factors. 

Motivated by large 03 (or 7) hints |2^, we show the 03-dependence of the branching 
fractions for /2 = (^2 = 5/2 in Fig. ^ The larger rates for larger 03 come from tree-penguin 
interference as in the K~7r+ case 112011 . 



B. B ^ VV^i PpK 

Unlike the 5^ K{TP)pTx^ case, the decay amplitude of pph^ with h = n or K 

contains both the current-produced {J') and transition (T) contributions: 



M{pph ) = J{pph ) + T{pph ) , 



(22) 



where 

Jipph" 



X (^pp 



{h-\{qb)v-A\B-) 



VubV*ga2 iuu)y_^ - VtiMtg 



as {uu + dd + ss)y_^ + 04 {qq)y_^ 



{uu + dd + ss)yj^^ ~'~ 2 '^^ {^uUU + Cddd + e^ss) 



+2a,Vt,V:^{h-\{qh)s-p\B-){pp\{qq),^p |0) 



V-A 



(23) 
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with q = d or s for h = ir or K, and, from Eq. (|T]), 

Gp 



T{pph ) 



V2 



{VubV*a, - VtbV*a^){h-\{qu)v-A\0){pP\{ub)v-A\B-) 



-2a(iVtbV^g 



ml 



ih-\{qu)v,A\Q){vP\{uh)v+A\B-) 



mb{mu + Trig 
ah{pp\ub\B') + Ph{pp\u'j5b\B') 



(24) 



where 



VubV*q ai - VtbV*g I 04 ± ae 



2 ml 



(25) 



mi,{mg + m„ 

It is of interest to compare the above equations with the famihar two-meson decay am- 



phtudes [11|. For the B ^ pp transition part, the analogous transitions are B^ 
(or the isospin-related ^ 7r+, p+). To single out this effect, we can search for two- meson 

7r+, p+ transition dominated 



decay modes dominated by such transitions. For the B^ 
modes, we have 5° — > n^n^{K^) and B'^ — > p^n^{K^) having decay amplitude propor- 
tional to and Pn{K), respectively. For the B^ 7r°, p° transitions, we can find ax, 
(5.„{^K) in B^ — i> Ti^K~ , p^TT~{K~) decay amplitudes, respectively. The vr^vr^ mode is different 
due to the cancellation of strong penguin in B^ and B^ tt^ transition parts as 
they are related by isospin. 



For the current-produced part, we can find similar terms in B 



and p^K~ decay amplitudes. However, we have additional terms. The matrix element of 
the isosinglet currents {pp\{uu + dd + ss)v,a\^) is non-vanishing, in contrast to the two- 
body K^TT^{p^) and 7r^7r'^(p'^) cases, where ir^lp^), as a member of an isotriplet, cannot be 
produced via the isosinglet current. As we will see, {pp\{uu+dd+ss)v,A\0) give non-negligible 
contributions to the ppK^{K^) modes. 



1. B ^ ppTT 

Although we have both current-produced and transition contributions to the ppTT~ mode, 
the former is expected to be small due to color-suppression of the 02 tree contribution and 
the smallness of CKM-suppressed penguin contributions, as one can see from Eq. (^31) . We 
show in Table 0(a) the contribution from the current-produced part J{pp'n:^)- As in the 
previous section, we are interested in the iVc = 3 case and list other cases for estimation of 
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TABLE V: (a) Current-produced Bj{B pp-K ) in units of 10 ^, and (b) strength of transition 
coefficients Ca,p,V5 giving rise to Br{B~ —>■ ppn") = 1.9 x 10~^. 

(a) Bj{ppTT~) for (ps = 54.8° (90°) in units of lO'^ 



Bv Ba 


Nc = 3 
Bv Ba 


Nc = oo 
Bv Ba 


G^j^ 0.03 (0.05) 0.03 (0.03) 
Fv,Dv,Sv 0.05 (OM) — 
Fa,Da,Sa — 0.09 (0.09) 


0.02 (0.03) 0.01 (0.01) 
0.00 (0.00) — 
— 0.02 (0.03) 


0.03 (0.03) 0.09 (0.12) 
0.15 (0.10) — 
— 0.29 (0.36) 


(b) \Cx\ values for cp^ = 


54.8° (90°) giving Bt -- 


= 1.9 X 10"^ 




2 Nc = 3 


Nc = oo 



\Ca\ (GeV^) 56.61 (63.52) 53.28 (59.86) 47.66 (53.66) 
\Cp\ (GeV^) 1233 (1356) 1160 (1278) 1038 (1146) 
ICysl (GeV^) 57.53 (57.53) 54.11 (54.26) 48.36 (48.77) 



non-factorizable effects. For = 3 the main contribution is from the strong penguin terms 
(04, ae), while the tree contribution is small due to the smallness of 02- For Nc = 2, 00, we 
have larger 02 and the main contributions come form the tree amplitude. With or without 
non-factorizable parts, the current-produced contribution is indeed much smaller than the 
experimental rate B {ppn-) = (1.9+^;° ± 0.3) x lO^^ [||. 

Since the current-produced part gives small contribution and the transition part T is gov- 
erned by tti, we expect the latter to give major contribution in the ppTc rate. The transition 
part T involves unknown B~ — > pp transition form factors -Fa.p.vs in the matrix elements 
{pp\{ub)s,p\B~) . We illustrate with three cases where only one form factors dominates. In 
each case we fit the coefficient Ca, Cp or Cys to the central value of the experimental mea- 
sured pp7r~ rate. Note that the matrix elements {pp\{ub)s,p\B~) in Eq. (p^ have nothing 
to do with the factorized meson h, hence the obtained coefficients Ca,p,v5 can be applied to 
the ppK~ mode (and for the ppK^ mode through Eq. (^8|)) as well. 

Table 0(b) shows the obtained values of these coefficients. It is interesting to observe 
that CAyb ~ {j^bI'^^ and Cp ~ [jnBl'if ■ Note that the effect of Fa and Fys in the pp'n 
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TABLE VI: B{B- pp-n-) in units of 10^^ for (/.a = 54.8° (90°) with h = d2 = z/2, S2 = S2 = 0. 







Nc 


= 2 


Nc 


= 3 


Nc 


= oo 






+\Cx\ 




+\Cx\ 


-\Cx\ 


+\Cx\ 


-\Cx\ 


X 


= A 


2.24 (2.42) 


1.84 (1.66) 


1.88 (2.07) 


1.96 (1.79) 


1.83 (2.01) 


2.85 (2.70) 


X 


= P 


1.99 (2.15) 


2.10 (1.92) 


1.86 (2.06) 


1.97 (1.80) 


2.29 (2.50) 


2.39 (2.20) 


X -- 


= V5 


2.62 (2.55) 


1.47 (1.53) 


1.99 (2.03) 


1.85 (1.82) 


1.41 (1.64) 


3.27 (3.06) 




2 2.5 3 3.5 4 4.5 5 2 2.5 3 3.5 4 4.5 5 

mpp (GeV) rripp ( GeV) 



FIG. 4: dB/dnipp spectrum for the B — > ppvr mode with Nc = 2> and (j)^ = 54.8°. Sohd, dash 
and dot-dash hues in (a) and (b) correspond to =F|C'y|) ^1^*^! ^"^^ T|C'p|, respectively. 

decay rate are similar. For </)3 = 54.8° (90°), we have \a.„\ ~ (0.9) 1/3,^1 for the color allowed 
tree-dominated part, leading to \Ca\ ~ (1-1) \Cv5\- Unlike the current-produced part, the 
transition part is not sensitive to Nc^ since a^, [3t, (composed of ai, 04 and ag) do not depend 
strongly on Nc as 02. 

By combining the transition contributions with the current-produced part for /2 = ^2 = 
5/2, S2 = S2 = case, we obtain the total branching fractions shown in Table [V^. We see 
that the rates are still similar to the experimental central value, justifying our procedure. 

In Fig. §we plot dB/dnipp for the B~ — * ppTr~ mode, for the Nc = 3 case with 03 = 54.8°, 
Ca = =F|Ca|, Cp = =f|Cp| and Cv5 = t\Cv5\, respectively. It is clear that the three cases 
give close to identical results. The threshold enhancement phenomena is evident, as has 
already been shown in the b'^ A(S°)p7r+ cases. However, we have a much faster l/rupp 
suppression here due to the behavior of the dominant transition form factors, while 

for the A(S°)p7r+ modes the form factors only behave as l/t"^ in the large t limit. It would 
be interesting to verify the faster l/rUpp fall off experimentally. 
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FIG. 5: (/)3-dependence of i? ppir branching fraction for /2 = ^2 = z/2, S2 = §2 = and 
Nc = 3. Solid, dash and dot-dash hnes for — |Ca|, — |Cp| and — ICysl, respectively. The plots are 
fixed to (ps = (a) 54.8°, (b) 90° values given in Table ^ 

In Fig. ^we illustrate the 03-dependence of the ppTT~ branching fractions, with 03 = 54.8° 
and 90° results fixed to those of Table |V|. Since this mode is dominated by the transition 
part, we expect similar 03 behavior as in tt+tt^, p+vr^ decay rates. The behavior of 

the rates from the Ca, Cp terms are similar to the tt+tt^ case |^0| rather than the 7r~7r° 
case. We have a similar term in the tt+tt" amplitude, while due to cancellation in the 
current-produced and the transition terms, there is no strong penguin in the 7r~7r° amplitude, 
resulting in small tree-penguin interference. The Cv5 case is similar to p+vr^ and does 

not show strong 03 dependence. As noted before, both features can be understood from the 
expression of a-,, and P-,, in Eq. (|25|). 



2. B ^ ppK 

For the B~ —>■ ppK~ decay, we have both a transition (T) part and now a more effective 
current-produced {J') part, as shown in Table VII| (a). For the vector part, the largest 



contributions come from (03 + a5){pp\{uu + dd + ss)v\0) and a4{pp\{ss)v\0) in Eq. (^31). In 
the first line of the table, we have {pp\ {ss)v,a,s,p\0) = and the contributions are mainly from 
the 03 + 05 term. For Nc = 3, 03 + 05 is small, resulting in a small By- In the second line, the 
04(^^1 (ss)v|0) term gives By ~ 0.3 x 10~^ and interferes differently with the previous term for 
different Nc as 03 + 05 changes sign. On the other hand, the (03 — a5){pp\{uu + dd+ ss)a\0) 
term dominates in the axial part. The dependence on of these contributions can be 
understood from the behavior of 03 — 05. 
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TABLE VII: (a) Current Bj{B~ ppR-) and (b) transition Bj{B- ppR-) in units of 10"^ . 



(a) Bj{ppK-) in units of lO'^ for (/-g = 54.8° (90°) 



Nc = 2 


Nc = 3 


Nc = oo 


Bv Ba 


Bv Ba 


Bv Ba 


G^j^ 0.02 (0.02) 0.10 (0.09) 


0.01 (0.01) 0.06 (0.06) 


0.33 (0.31) 0.01 (0.02) 


Fv,Dv,Sv 0.65 {Om) — 


0.26 (0.26) ~ 


0.02 (0.03) — 


Fa,Da,Sa — 0.23(0.19) 


— 0.12 (0.12) 


— 0.02 (0.05) 



(b) BrippK-) in units of 10"^ for 03 = 54.8° (90°) 

Nc = 2 iVc = 3 Nc = oo 

Ca 2.08 (3.42) 2.12 (3.48) 2.18 (3.58) 

Cp 1.83 (2.86) 1.86 (2.91) 1.90 (2.99) 

Cv5 0.26 (0.20) 0.23 (0.19) 0.18 (0.18) 



We now turn to the transition (T) part. For 03 ~ 54.8° (90°), we have \aK\ ~ |«7r| 
{Wk\ ^ W-k\) and hence Br{ppK~) ~ (or >) Br{j>pT^~) from the Ca,p contributions. On 
the other hand, 04 and ag in /5x are partially canceled and the tree contribution is CKM 
suppressed, resulting in -C Hence, the {pp\uh\B^) contribution containing 

Ca,p is much larger than the {pp\u'^^h\B~) case coming from Cys, as can be seen from 
Table |VIl| (b). Note further that the 03 = 90° case from Ca,p gives larger rates, as should 
be expected from the analogous K^ti~ mode ||20| . 

We combine the current-produced contribution with the transition part and give the total 
branching fractions in Table [VII I| . The results prefer the 03 = 90° case. Numbers shown in 
the first two lines of Table VIII are close to the the experimental rate B{B^ ppK^) = 



(4.3+^:^ ± 0.5) X 10-6 [|. We see that, for the 03 = 90°, = 3 and Ca = -\Ca\ case we 
have B = 4.26 x 10"^, which is closest to the central value of the experimental rate. The 
value only changes by 10% as we modify Nc to 2 or 00. 

We plot dB/drripp in Fig. ^for the 03 = 90°, Nc = 3 case of Table VIII . The curves from 



Ca, Cp terms are close to data points taken from Ref. |]T|. The curve from Cv5 term is too 
low as expected from the smallness of Pk- For the first two cases, one can see that, except 
for a possible bump at nipp ~ 2.2-2.4 GeV, the behavior of the decay spectrum including 
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TABLE VIII: B{B- ppR-) in units of lO^^ for 03 = 54.8° (90°), with /s = da = z/2, 

S2 = S2 = 0. 





N, = 2 


Nc = 3 


Nc = oo 




+\Cx\ -\Cx\ 


+\Cx\ -\Cx\ 


+\Cx\ -\Cx\ 


X 


= A 2.56 (3.85) 3.36 (4.75) 


2.18 (3.45) 2.82 (4.26) 


2.06 (3.36) 2.38 (3.96) 


X 


= P 2.72 (3.75) 2.70 (3.73) 


2.25 (3.30) 2.23 (3.28) 


1.95 (3.08) 1.94 (3.07) 


X -- 


= V5 0.83 (1.08) 1.45 (1.08) 


0.42 (0.57) 0.80 (0.56) 


0.16 (0.16) 0.28 (0.36) 




FIG. 6: dB/drripp of the ppK~ mode with (ps = 90° and Nc = 3. The sohd, dashed and the 
dot-dashed hnes in (a) and (b) stand for the Ca = ^fIC'aI, Cp = =f|Cp| and the Cv5 = T\Cv5\ 
case, respectively. 



threshold enhancement can be explained naturally. Comparing with Fig. ^, we note that 
the suppression for large rripp is milder than the ppn" case. This is due to the presence of 
the current-produced part which has less suppressed form factors (1/t^), which dominates 
over the transition part in the large t region. 

We give in Fig. |^the 03-dependence of B{ppK~). The Ca,p cases are similar to the K~tx^ 
mode while the Cys case is similar to the case as discussed before. We note that 

the experimental indication that the ppK~ rate is larger than the ppTr~ rate seems to favor 
larger 03 values such as 90° case, analogous to B to two meson decay situation . 
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FIG. 7: 03-dependence of B{B ppK ). Solid, dash and dot-dash Unes stand for the Ca 
— \Ca\, Cp = —\Cp\ and Cv5 = — |Cv5| cases, respectively. 



3. B 



■0 



ppK^ 



For the B ppK^ decay, we have M.ippK'^) = JippK^) + T[ppK^). Since 

{W\{sh)v-A\B^) = {K-\{sb)v-A\B-) by isospin, we have 



J{ppK'>) = J{ppK~ 



On the other hand, 



T{ppK^) = i ^f-RO ruh 



(26) 



(27) 



with 



a 



mbims + md) 



(28) 



We show in Table |^ the separate current-produced and transition contributions to the 
ppK^ decay rate. As explained in the above, the current-produced part is identical to the 
ppK~ case, except for the difference in Tp- and r^o. For the transition part, the transition 
form factors are related to the ppK~ case through Eq. (plsl) . We concentrate on the {db)s 
form factors instead of {db)p, since la^ol ^ |/%o|. We have J-'a = I-I-^a + 0.9Fv5 and 
jFp = —Fp/A for the {db)s form factors. Since = —1/4, the contribution form 

the Cp term is very small. The ratio of the Ca and Cy^ contributions can be understood 
as well. For 03 = 54.8°, the Fa dominated case is larger than the Fy^ dominated case by 
a factor of 11/9 in amplitude, giving a rate enhancement ~ (Il/9)^|CA/Cy5p ~ (11/9)^ ~ 
1.5. For 03 = 90°, we have a further 10% growth in amplitude due to \Ca\ ~ l-ljCysl 
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TABLE IX: (a) Current Bj{B ppK°) and (b) transition Bt(B^ ppK°) in units of lO"^. 
(a) Bj{pplf) in units of 10"^ for ^3 = 54.8° (90°). 





Nc = 2 


Nc = 3 


Nc = 00 




Bv Ba 


Bv Ba 


Bv Ba 


G^" 0.02 (0.02) 0.10 (0.08) 


0.01 (0.01) 0.06 (0.05) 


0.31 (0.29) 0.01 (0.02) 


Fv,Dv,Sv 0.61 (0.65) — 


0.25 (0.24) — 


0.02 (0.03) — 


Fa, Da, Sa 


— 0.21 (0.18) 


— 0.12 (0.11) 


— 0.02 (0.05) 


(b) Br{ppK^) in units of 10-^ for (p^i = 54.8° (90°). 




Nc = 2 


A^c = 3 


iVc = 00 


Ca 


2.74 (3.37) 


2.78 (3.43) 


2.84 (3.52) 


Cp 


0.12 (0.14) 


0.12 (0.15) 


0.13 (0.15) 


Cv5 


1.91 (1.87) 


1.93 (1.89) 


1.96 (1.95) 



(as shown in Table 0(b)), and the rate enhancement becomes ~ 1.8. It is interesting to 
compare the transition contributions to those in the ppK~ mode: for the Ca dominated 
case, BrippK^) ~ Br{ppK~); for the Cp dominated case, BrijppK^) ^ Br{ppK~); for the 
Cv5 dominated case, BrippK^) ^ Br{ppK^). 

We give in Table ^ the full branching fraction by combining the current-produced and 
transition parts in amplitude. For the = 3, 03 = 54.8° [90°] case, we have B{ppK^) = 
(0.5-3.6) X 10~^ [(0.5-4.3) x 10"^]. It could be close to or smaller than the ppK~ rate. 

In Fig. I we plot dB/dnipp for the 5^ ppK'^ mode. The decay spectrum for the Ca 
and Cys cases are similar to the Ca, Cp cases in the ppK~ mode. They exhibit threshold 
enhancement and a slower fall off for large t compared to the ppTi' case. For the Ca case, 
where the rate is not far from ppK~ , the decay spectrum could be checked soon. 

In Fig. H we show the 03-dependence of B(j>pK^). As shown in Eq. (|^), the transition 
part does not have tree-penguin interference and hence the mild 03 dependence is from the 
sub-dominant current-produced term. 
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TABLE X: B{B- ppK^) in units of 10"^ for ^3 = 54.8° (90°) and h = d2 = z/2, S2 = §2 = 0. 





Nr = 2 


Nr = 3 


A^'c = 00 




+\Cx\ -\Cx\ 


+\Cx\ -\Cx\ 


+\Cx\ -\Cx\ 


X 


= A 2.91 (3.53) 4.21 (4.86) 


2.68 (3.28) 3.60 (4.28) 


2.71 (3.32) 3.05 (3.88) 


X 


= P 0.95 (0.97) 0.94 (0.96) 


0.48 (0.50) 0.48 (0.49) 


0.17 (0.23) 0.16 (0.22) 


X -- 


= V5 2.13 (2.14) 3.33 (3.25) 


1.88 (1.85) 2.70 (2.64) 


1.86 (1.82) 2.14 (2.24) 




uipp (GeV) mpp (GeV) 

FIG. 8: dB/drUpp for the ^ — > ppK^ mode with /2 = ^2 = z/2, S2 = §2- Solid, dash and dot-dash 
lines in (a) and (b) are for Ca = =F|C'a|, Cp = =f|C'p| and Cvb = T\Cv5\ cases, respectively. 
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FIG. 9: 03-dependence of B{B'' —>■ ppK^) for 72 = ^2 = z/2, S2 = §2 = 0. Sohd, dash and 
dot-dash lines are for Ca = —\Ca\, Cp = —\Cp\ and Cv5 = —\Cv5\ cases, respectively. 

C. Comparison with Other Works 

Before we end this section, we compare our work with some others. There are approaches 



that use pole models to evaluate decay matrix elements. For example, Ref. [21 1 use K* pole 
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for Kp production in B rj'Ap decay, while Ref. |^ use pole models in two-body and 
three-body baryonic B decay. We focus on the comparison with the Ref. as we have 
some subjects in common. 

In Ref. Cheng and Yang use a factorization approach in the current-produced {J') 
amplitude. Their approach is similar to ours (up to some technical differences), hence 
they obtain results similar to ours in current production dominated modes, such as 5° — >■ 
Apn^, S°p7r+. However, there is considerable difference in the transition (T) part. We 
factorize the amplitude into a current-produced meson and a B to baryonic pair transition 
amplitude. In their approach, they use a simple pole model to evaluate this part. For 
example, in B^ ppK^ decay, they have a strong process B- ^ {aI*\ Sj(*)}p, followed by 
a weak {A^*'', S^*-*'*} pK~ decay. From their modeling of the strong coupling gA,,-^B-p = 
Sv^Oso^B-o) they can give ppK~ rate that is close to experimental result by using monopole 

b ^ 

dependence of g\^^B-p- The nipp spectrum given in Ref. shows a peak around t ~ 
6 GeV^ (or rupp ~ 2.5 GeV), while we have a sharper peak in lower nipp (around 2 GeV). The 
difference is due to the behavior in our transition part from QCD counting rule, while 
they have ~ from the pole model. On the other hand, one expects peaking behavior 
towards large mpx- due to Ab pole in their approach, while in this work we do not expect 
any structure (since K~ and p are factorized) in the nipK- spectrum. In turn, they expect 
ppK^ rate ~ 10~^ due to the absence of A;, pole, while we expect a rate that could be as 
large as ppK~ , although 10"'' is also possible. It is up to experiment to check the rripx- 
spectrum and the ppK'^ rate. 

We recall that in Ref. |^ we also tried a VMD (dispersion analysis) approach |]15|, ^ in 
the current production dominated B^ D*"pn decay. In this approach, the strong coupling 
for each pole is fixed at the pole mass, hence each pole gives a monopole contribution to 
the total form factors. One needs to have more than one pole with cancellations in order to 
reproduce the correct QCD counting rule, which is for current-produced form factors 
||15| , p^ . We likely would have the same situation here, that more than one pole is needed to 
reproduced the large t behavior. If we take a multi-pole approach, the baryonic transition 
form factor can be expressed as B ^ transitions with Mj as one of the mesons, followed 
by a strong process Mj BB (similar to Ref. |21|). Summing over i, the QCD counting rule 
should be taken as a constraint. Instead of doing so, in part because of lack of independent 
data, we use a simplified transition form factor motivated from the QCD counting rule 
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directly in this work, and wait for experimental data, such as semi-leptonic B BB'lu and 



semi-inclusive B BB'X (similar to i? ^ nX |2^), to improve our understanding. One 
may resort to the multi-pole approach once these measurements become available. 



V. DISCUSSION AND CONCLUSION 



In this work we use factorization approach to study charmless three-body baryonic B 
decays. We apply SU(3) relations and QCD counting rules on baryon form factors. We 
identify two mechanisms of baryon pair production, namely current-produced and transition. 
The Il^p7T~^ and ApvT"*" modes arise solely from the current-produced part, with rates of order 
10^^ — 10^^. The ppTT^ , ppK^ and ppK^ modes are dominated by transition contributions, 
while the current -produced contributions in the last two cases are significant. 

Due to the absence of Sv,a in the current-produced amplitude, and the complete absence 
of the transition amplitude, 'E^pn'^ and Apvr^ modes are the simplest in this work. However, 
they are sensitive to how we treat the chiral limit of the pseudoscalar term (which has an ag 
coefficient). On the other hand, we neglect F2 contribution in the vector part. It remains to 
be checked whether this is a good approximation or not. It may in turn give us information 
on F2, or equivalently Ge, from these measurements. In particular, the vector current form 
factors in ApvT"*" are only related to the proton Gm,e from SU(3) symmetry (as one can see 
from Eq. ( ^I]) and Table p| ), and we may obtain information of from this mode. 

The 03 dependence of Apyr"*", E'^pvr"'" rates are similar to K^tc^, K*^tt^ modes. Since 
transition contributions dominate in the pph modes, the 03 dependence of ppn~{K^) rate 
is similar to that in the two-body B^ — > -n-^n^ (K^) or p^n^i^K^) decays, while the 03 
dependence of ppK^ rate is mild. 

Under factorization, the ppT{~ and the ppK~ modes have the same baryonic transition 
form factors. Since the ppTt" mode is dominated by baryonic transition contribution, we 
use it to fit for the transition form factor parameters and apply them to the ppK~ case. To 
keep ppTT~ rate around the experimental central value but allowing ppK~ rate to be larger, 
data seems to favor a larger 03. The ppK^ rate can be similar to or much smaller than the 
ppK~ rate. We do not consider modes involving vector mesons, such as ppK*, since they 
will involve further unknown form factors. 

It is interesting that we can reproduce ppK~ decay spectrum based on QCD counting 
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rules, indicating that the latter is a rather robust theoretical tool. From QCD counting rules, 
the current-produced baryonic form factors behave like while the transition baryonic 

form factors behave like in large t (= rn^^,) limit. The ppTi~ decay is dominated by 
the transition part and hence shows a faster damping behavior for large t. In contrast, the 
Ap7r+, S°p7r+, ppK^ and ppK^ modes contain current-produced part, and show a slower 
damping behavior for large t. We expect an even slower damping behavior (form factors 
~ 1/t) in three-body mesonic B decay. These can be checked soon, especially by comparing 
the ppK' and the ppix^ spectra. 

We note that there is a possible bump around rupp = 2.3 GeV in the ppK~ decay spec- 
trum. It is interesting that the position is close to the mass of a glueball candidate ^(2230), 
also known as /j(2220), with m = 2231 ± 3.5 MeV and T = 23^7 MeV [||]. Combining 
B{J/ip 70 > 2.5 X 10"3 [H] and B{J/ij ^ 7O ^ pp) = (1.51^;^ ± 0.5) x 10"^ |5 



we have pp) < 6 x 10 ^. If the rate in the whole 2.2 < Mpp < 2.4 GeV bin is due 

to this resonance, we would have B{B^ K^^)B{^ pp) ~ 1.24 x 10~^. We thus get 
B{B- K-^) > 2xlO-^orafewtimestherateof;B(5- r]'K-) = (6.5±1.7) x 10^^ [|4 



Since both ^ and 77' are glue rich hadrons, while b ^ s decays provide glue rich environ- 
ment this may be of great interest. The underlying dynamics could be g* g^, which 
is analogous to g* — > gr]' for B ^ rj' + Xg decay . One should also search in three-body 
mesonic decay modes. However, as noted, a slower fall off for non-resonance part, together 
with interference with possible nearby resonances, may produce physical background. But 
B — > ppK^ decay could be a rather clean mode to search for ^ p8| . 
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APPENDIX A: SOME USEFUL FORMULAS 



In general, for a three-body decay B /iBB , the amphtude can always be written in 
the following form: 



M (B ^ hBB' ) = % <{^m(pb)A^(Pb') + 'Sw(Pb)M5^(Pb'' 



V2 



(Al) 



whose absolute square is given by 



2,2 2 2\/2 2 2 

+ m2 - - "^23 "^23 - "^2 " "^3 



+ 2Re{AV* 



nil (^^23 ~ ^2 ~ ~^ ""^2 (^TTT'^B + ^2 ^ 



m 



23 



12 



(mi + 



i3, P ^ C, ma 



-ma) 



(A2) 



where the summation is over all spins, and we have adopted the convention that the relatively 
positive baryon is assigned as particle 1, the other baryon is assigned as particle 2, the 
meson h is always assigned as particle 3. One can see from the above that only Iie{AT>*) 
and Re(i3C*) appear as interference terms upon squaring. Given these formulas, the task is 
now reduced to obtain the A-V terms for an amplitude of interest. 

It is straightforward to obtain the decay rate F from the integration of 

1 1 



dT 



{2 7rf 32 m 



(A3) 



B 
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